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چکیده 


تغییر اقلیم اثرات قابل توجهی بر چرخه آبی دارد» به گونه‌ای که می‌تواند تهدیدی بزرگ برای سامانه‌های آبی در جهان قلمداد شود. هدف از این 
مطالعه بررسی اثر احتمالی تغییر متفیرهای اقلیمی بر رواناب و تبخیر و تعرق حوزه ابخیز مهرگرد است. برای این منظور با استفاده از مدل 91951715.2 
و سناریوی ول ۵ ۱۱۵۱ خروجی مدل تعییر اقلیم ٩0۳22‏ ریزمقیاس‌نمایی و برای دوره ۷ :۲۰۳۲۰ شبیه‌سازی شد. از داده‌های روزانه حداقل و 


حداکثر دما در ایستگاه بروجن و بارش ایستگاه تنگ زردآلو برای دوره پایه ۱۹۸۴-۲۰۰۵ به عنوان ورودی‌های مدل اقلیمی استفاده شد. هم‌چنین برای 
شبیه‌سازی شرایط هیدرولوژیکی حوزه آبخیز از مدل ۸2012 5۷۷ استفاده شد. عملیات واسنجی برای دوره ۲۰۰۴-۲۰۱۲ و اعتبارسنجی برای دوره 
۲۰۱۳-۲۰۱۶ با استفاده از الگوریتم 517۳1-2 انجام شد. نتایج حاصل از ارزیابی مدل 5195[15.2 براساس معیارهای آماری 5 12 و 11512 
نشان داد که مدل در شبیه‌سازی متغیر دما نسبت به بارش از دقت بالاتری برخوردار بوده است. بر این اساس برای آینده کاهش ۵۳/۴۸ درصدی بارش و 
افزایش ۰/۸۴ و ۳/۹۹ درجه سانتی‌گراد به‌ترتیب برای حداقل و حداکثر دما پیش‌بینی شد. ضریب نش-ساتکلیف و ضریب تعیین برای رواناب در مرحله 
واسنجی به‌ترتیب ۰/۶۹ و ۰/۷۳ و برای مرحله اعتبارسنجی به‌ترتیب ۰/۵۸ و ۰/۷۱ به دست آمدند. هم‌چنین تغبيراقليم باعث کاهش ۲۳/۸۲ درصدی 
رواناب» کاهش ۲۶/۰۳ درصدی تبخیر و تعرق واقعی و افزایش ۱۰/۲۰ درصدی تبخیر و تعرق پتانسیل خواهد شد. 


واژه‌های کلیدی: بیلان آب. تغییر اقلیم» ریزمقیاس‌نمایی» سناریوهای اقلیمی» مدل 5۱۷۸۲ 


مقدمه 

اقلیم کره زمین از زمان پیدایش جو بر روی این سیاره تغییر 
کرده» اما سرعت آن هیچ گاه به اندازه عصر حاضر نبوده است. از نظر 
پژوهش گران» اين مقدار شدت در تغییرات اقلیمی به‌دلیل افزایش 
گازهای گلخانه‌ای در اتمسفر است که غالبا منشأً انسانی دارند (۱۶). 
بر اساس گزارش ارزیابی پنجم (۸185) هیثت بین‌الدول تغییر آب و 
هوا (1۳06» دمای جهانی طی سال‌های ۱۸۸۰ تا ۲۰۱۲ به میزان 
۵ «درجه سانتی‌گراد افزایش يافته است و در آینده تا سال ۲۱۰۰ به 
میزان ۰/۳ تا ۴/۸ درجه سانتی‌گراد افزایش خواهد یافت (۱۷). تغییرات 
آب و هوایی از طریق تغییر در زمان» مقداره شکل بارش, نرخ تبخیر و 


۱ و ۲- به‌ترتیب دانش آموخته کارشناسی ارشد مهندسی آبخیزداری و استاه گروه 
مرتع و آبخیزداری» دانشکده منابع طبیعی دانشگاه یزد 
(*+- نویسنده مسئول: ۴ 210151 :اتفصص۴) 
۲- دکتری ژئومورفولوژی مدیر دفتر برنامه‌ریزی و بررسی‌های اقتصادی شرکت آب 
منطقه‌ای اصفیان 
09 :]10 
صمت مافصتل مه امصوط آم)ممصصه0۷عع]18 -4 


انتقال و رطوبت خاک» رواناب بر چرخه هیدرولوژیکی تأثیر می‌گذارد 
(۲۳). استفاده از مدل‌های هیدرولوژیکی همراه با مدل‌های گردش 
عمومی (015) و یا مدل‌های آب و هوایی منطقه‌ای (615) 
رایج‌ترین روش برای ارزیابی تأثیرات تغییرات آب و هوایی بر جریان 
است (۲۵ و ۸) در مطالعات هیدرولوژی باید توجه داشت که معمولا 
مساحت سلول‌های محاسباتی مدل‌های ۸000۲" بزرگ‌تر از 
وسعت منطقه مطالعاتی است بنابراین برای استفاده از نتایج این 
مدل‌ها در سطح حوضه باید خروجی‌های مدل ریزمقیاس نمایی شوند 
(۱۴). روش‌های مختلفی بدین منظور وجود دارد که روش‌های آماری 
از روش‌های مرسوم می‌باشد و نرم‌افزارهایی نیز در این رابطه توسعه 
یافته است که از جمله آن‌ها می‌توان مدل 5۳5[۷ ۸ را نام برد (۱۵). 
شناخت و درک مفهوم بیلان آب یکی از پیش نیازهای ضروری در 
راستای مدیریت پایدار منابع آب در سطح حوزه‌های آبخیز است. 
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مهمترین عناصر بیلان آب از دیدگاه مدیریت منابع آب شامل بارش» 
تبخیر و تعرق» رواناب سطحی. جربان زیرزمینی و جریان جانبی 
می‌باشند. بنابراین» می‌توان اجزاء منابع آب موجود در یک سیستم را 
در دوره‌های زمانی مختلف مورد مقایسه قرار داد و درجه تأثیر هر 
یک را در تغییرات بیلان آبی سیستم‌های هیدرولوژیکی شناسایی کرد 
(۲۱). امروزه برای مطالعه و برنامه‌ریزی پایدار و موّثر در مدیریت 
جامع حوزه‌های آبخیز از مدل‌های هیدرولوژیکی استفاده می‌شود. 
ابزار ارزیابی آب و خاک (5۷۷۸۲) نمونه‌ای از مدل‌های 
هیدرولوژیک با مبنای فیزیکی است که بر اساس ویژگی‌های حوزه 
آبخیز و شرایط اقلیمی آن به شبیه‌سازی در مقیاس بزرگ و هم‌چنین 
پایش فرایندهای مرتبط با چرخه آب می‌پردازد (۲۹). با توجه به 
اهمیت بحث اثر تغییر اقلیم بر منابع آب در این زمینه تحقیقات زیادی 
انجام شده است. زو و همکاران (۳۸) در مطالعه‌ای پاسخ هیدرولوژیکی 
به تغییرات آب و هوایی در حوضه رودخانه ۷۷۵ چین را با استفاده از 
مدل 5۷۷۸۲ و ٩25۷۲‏ مورد بررسی قرار دادند. نتایج نشان داد که 
میانگین رواناب سالانه در دوره ۲۰۴۶-۲۰۶۵ ۲۰۸۱-۲۱۰۰ به‌ترتیب ۱۳/۴ 
درصد و ۴۵ درصد نسبت به دوره پایه ۱۹۶۱-۲۰۰۸ افزایش می‌بابد. ژانگ 
و همکاران (۳۷) در مطلعه‌ای با استفاده از دو مدل 50/۸۲ و 657 
جریان‌های رودخانه 260 در چین را تحت تأثیر تغییرات اقلیمی 
ارزیابی کردند. در ادامه از 00۷های مختلف (1300-05(1.1 
2 و 0۲۳5[۷۲1-۷۲() تحت سه سناریو (0۳5) استفاده 
کردند. بر این اساس درجه‌حرارت به‌طور مداوم در منطقه افزایش 
می‌یاب با این حال بارش آینده بسیار پیچیده و نامطمئن است. ملک 
و آبگاز (۲۲) در مطالعه‌ای اثر تغییر اقلیم بر شرایط هیدرولوژیکی 
حوضه رودخانه 0۷۳272 را با استفاده از مدل هیدرولوژیکی 5۷۷۸۲ 
و ٩۳5[۷‏ مورد بررسی قرار دادند. نتایج نشان‌دهنده آن است که 
میانگین جریان سالانه در آینده به میزان ۱۷/۸ درصد افزایش می‌بابد. 
آبراهام و همکاران (۲) در مطالعه‌ای شرایط رواناب آینده تحت تأثیر 
تغییر اقلیم را توسط دو مدل اقلیمی 01۷۲۳5 در دریاچه 717۷2۷7 دره 
ریفت اتیوپی مورد بررسی قرار دادند. نتایج نشان داد که حداقل و 
حداکثر دما تحت سناریوهای 0۳8.5 و 130۳4.5 افزایش می‌یابد و 
در مقابل بارش به شدت کاهش می‌یابد. قرمزچشمه و همکاران (۱۰) 
در پژوهشی تغییرات جریان تحت تأثیر اقلیم در حوضه رودخانه کن را 
با استفاده از مدل هیدرولوژیکی 5۷۷۸ و مدل 
ریزمقیاس‌نمایی ٩195‏ مورد بررسی قرار دادند. در نهایت برای همه 
ایستگاه‌های مطالعاتی افزايش دمای بيشینه و کمینه به ترتیب حدود 
۳ و ۰/۸ درجه سانتی‌گراد و کاهش بارش پیش‌بینی شد. نتایج 
حاصل از شبیه‌سازی رواناب» کاهش آن را نسبت به دوره پایه نشان 
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داد. اکبری و همکاران (۴) در پژوهشی اثر تغییر اقلیم بر رواناب 
شضه روا ای زا یمه فل ‏ و ستاوها 
گزارش پنجم 1۳0 مورد بررسی قرار دادند. نتایج نشان‌دهنده 
کاهش رواناب در سال‌های بررسی شده است. گودرزی و همکاران 
(۱۳) در مطالعه‌ای عدم قطعیت روش‌های ریزمقیاس‌نمایی در تخمین 
جربان آتی حوضه آبریز قره‌سو تحت اثر تغییر اقلیم را بررسی کردند. 
نتایج حاکی از برتری مدل‌های 020۳5۷۲2 و ۲1200۳03 برای هر دو 
متفیر بارش و دما بود. در نهایت نتایج با هدف بررسی تغییرات رواناب 
تحت تآثیر تغییر اقلیم به مدل هیدرولوژیکی معرفی شد. به طور کل 
اکثر مدل‌ها افزایش رواناب در فصل زمستان و کاهش آن در بقیه 
فصول نسبت به دوره پایه را پیش‌بینی می‌کنند. نظر به این که حوزه 
آبخیز مهرگرد در سال‌های گذشته در وضعیت مطلوبی به لحاظ منابع 
آبی قرار داشته است اما متأسفانه در چند سال اخیر به‌واسطه تغییر 
شرایط آب و هوایی» وضعیت هیدرولوژیکی آن دستخوش تغییراتی 
شده است به‌نحوی که ساکنین حوزه آبخیز مذکور را به لحاظ 
دسترسی به منابع آب با مشکل رو به رو کرده است. بنابراین با توجه 
به اهمیت موضوع تغبیر اقلیم و تأثیر آن بر اجزاء مختلف کره زمین 
به خصوص چرخه هیدرولوژیکی» این پژوهش در نظر دارد تا با استفاده 
از مدل هیدرولوژیکی 5۷۷۸۲ و مدل اقلیمی 020۳61۷2 که یکی 
از مدل‌های تغییر اقليم جفت شده 0۷۲1۳5 است تحت سناریو 
5 اثر تغییر متفیرهای اقلیمی را بر رواناب و تبخیر و تعرق 
واقعی و پتانسیل حوزه آبخیز مهرگرد در آینده مورد بررسی قرار دهد. 
برای این منظور از مدل آماری 5795۷5.2 برای ریزهقیاس‌نمایی 
داده‌های اقلیمی بزرگ مقیاس پایگاه ۳0۳۲۳" استفاده شد. ضمن 
اینکه مدل 5151 تاکنون بیش‌تر برای ریز مقیاس‌نمایی دما و 
بارش در مناطق مرطوب به‌کار رفته است و در نواحی خشک و نیمه 
خشک کمتر مورد ارزیابی فرار گرفته است. لازم به ذکر است که در 
این حوزه آبخیز تاکنون هیچ‌گونه مطالعه‌ای از حیث بررسی اثر تغییر 
اقلیم بر مولفه‌های بیلان آب صورت نگرفته است لذا اين پژوهش 
می‌تواند زمینه‌ای را برای مطالعات گسترده‌تر در راستای مدیریت 
هرچه بهتر منابع آب در حوزه آبخیز مذکور و حوزه‌های آبخیز با 
شرایط مشابه فراهم کند. 


مواد و روش‌ها 
منطقه مورد مطالعه 


حوزه آبخیز مهرگرد به طول جغرافیایی ۵۱ درجه و ۴۷ دقیقه و 
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۱ درجه و ۵۰ دقيقه و ۵۷ انیه تا ۲۱ درجه و ۲۳ دقیقه و ۱٩‏ ثانیه 
شمالی به‌عنوان بخشی از حوضه آبریز کارون بزرگ با وسعت 
۹ کیلومتر مربع» در شمالی‌ترین بخش شهرستان سمیرم 
استان اصفهان واقع شده است (شکل ۱). کم‌ترین ارتفاع محدوده 
۵ متر مربوط به بخش‌های خروجی حوزه آبخیز و بیش‌ترین ارتفاع 
محدوده ۲۶۸۸ متر از سطح دریا مربوط به قسمت‌های شمالی می‌باشد. 
اقلیم حوزه آبخیز براساس اقلیم نمای آمبرژه نیمه خشک سرد و بر 
اساس روش دومارتن نیمه خشک می‌باشد. بارندگی متوسط دشت در 
این محدوده ۲۳۹ میلی‌متر است. هم‌چنین متوسط بارندگی در 
ارتفاعات ۲۶۸ میلی‌متر در سال براورد شده است (۲۴). 
شبیه‌سازی شرابط اقلیمی آینده با استفاده از مدل 
٩) (2‏ 

برای شبیه‌سازی شرایط اقلیمی آینده حوزه آبخیز مهرگرد از مدل 
2 استفاده شد. این مدل از جمله مدل‌های گردش عمومی 
جو سری 11۳5 می‌باشد و چهارمین نسل از مدل‌های آب و هوایی 
است که توسط مرکز مدل‌سازی و تحلیل آب و هوای کانادا زیر نظر 
سازمان محیط زیست این کشور توسعه یافته است. در اين مدل کل 
سطح زمین به‌صورت ۶۴*۱۲۸ سلول» شبکه‌بندی شده است (۶). 
هم‌چنین در پژوهش حاضر برای ریزمقیاس‌نمایی از مدل ٩151‏ 
انتفاده اد این مدل-یکی از مدل‌های. ریمقیاس‌نماین تابع اتقالی 
آماری به شمار می‌آید (۲۶). با توجه به‌اینکه در چند سال اخیر 


اثر تغییرات اقلیمی بر رواناب و تبخیر و تعرق حوزه آبخیز مهرگرد ۲۲۷ 


به‌واسطه تغییر شرایط آب و هوایی» وضعیت هیدرولوژیکی حوزه آبخیز 
مهرگرد دستخوش تغییراتی شده است مدل‌سازی میزان و نحوه 
تغییرات متفیرهای اقلیمی حوضه برای دوره آتی از اهمیت بالایی 
برخوردار است. بنابراین با توجه به هدف پژوهش, بازه زمانی ۲۰۳۰- 
۷ به عنوان آینده نزدیک مبنای بررسی اثر تغییر اقلیم بر رواناب و 
تبخیر و تعرق واقعی و پتانسیل حوضه قرار گرفت. همچنین به‌منظور 
بررسی اثر این تغییرات بر پارامترهای هیدرولوژیکی از داده‌های روزانه 
حداقل و حداکثر دما ایستگاه بروجن و بارش ایستگاه تنگ زرداآلو 
برای دوره پایه ۱۹۸۴-۲۰۰۵ به‌عنوان ورودی‌های مدل اقلیمی 
استفاده شد. بر این اساس با استفاده از مدل 8126315.2 : خروجی 
مدل گردش عمومی 020۳5[۷]2) به‌وسیله داده‌های مشاهداتی 
ریزمقیاس‌نمایی و با در نظر گرفتن سناریوی 1308.5 برای آینده 
شبیه‌سازی شد. سپس ۲۶ متفیر اقلیمی بزرگ مقیاس از پایگاه داده 
۳ در قالب یک فایل فشرده شده از طریق لینک 

وم8260100-2 02/102 :0282012 .و۲۱۵ متاقتاه 
۵ دانلود و به عنوان داده‌های ورودی به مدل 
اقلیمی معرفی شد. پس از انجام تکرار فرایند آنالیز از بین ۲۶ متغیر 
بزرگ مقیاس جو (پیش‌بینی‌کننده‌ها» متغیرهایی که بیشترین 
همبستگی را با بارش و دمای حوزه آبخیز (پیش‌بینی شونده‌ها) داشتند 
انتخاب شدند. برای اطمینان از صحت نتایج به‌دست آمده عملیات 
واسنجی و اعتبارسنجی انجام شد. 


ات سه نی 


حوزه آبخیز مهرگرد 
ار و ۱ 


وتو ۱ 
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سد مخزنی 


شکل ۱- موقعیت حوزه آبخیز مهرگرد در استان اصفهان و ایران 
و 200 0۳۵۲1۱6۵ هدهع از 0مصاوه)۲۵ ۱۷۲۵۵۰۵۵0 م6 ۵۶ حمتاهمنا -1 عساع1ز۲ 
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هم‌چنین برای ارزیابی کارایی مدل از یک سری معیارهای آماری 
مانند ۲5 ۲ ۳82 و ۳۱۷65۲" استفاده شد که در روابط ۱ تا ۳ به 


ترتیب ارائه شده است (۲۶). 


)۱ - 1 - ولا 
۳( - :00( 
1:21 
6۵-۵0-9 ۲ 
21<[ 2 
۳ 7 و۳9 


0 و 
۲-۳ 0-0 
ی( 0-0 
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| <- ر]۷]9 
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0 ۲ 
در روابط ۱ تا ۲ " مقادیر مشاهده شده " مقادیر شبیه‌سازی 
میانگیی مقافیه امه کم نگ قامن شا 
شده و 1 تعداد سال‌های آماری می‌باشد. 


شبیه‌سازی شرایط هیدرولوزیکی با استفاده از مدل 
52۸ 

مدل ٩۱۷۸۲‏ به لحاظ مقیاس فضایی نیمه توزیعی و از نظر 
زمانی مدلی پیوسته است. در مدل 51۷۸ ابتدا از روی مدل رقومی 
ارتفاعی (05۱8» حوزه آبخیز اصلی به تعدادی زیرحوضه تقسیم 
می‌شوده سپس بر مبنای نقشه‌های خاک و کاربری اراضی و شیب 
زیرحوضه‌ها نیز به واحدهای کوچکتری تقسیم می‌شونده که به 
هرکدام از اين واحدها یک واحد پاسخ هیدرولوژیکی (01181) 
می‌گویند. این واحدها باید تا حد امکان از نظر هیدرولوژیکی مشابه 
باشند. در واقع باید در هر واحد هیدرولوژیکی ویژگی‌های خاک» 
توپوگرافی» پوشش و کاربری اراضی تغییرات چشمگیری نداشته 


۱۳۱ 
عمتاحصتصهعاه 1۲ 0۶ اصمام006111) -2 

3- 1500] ۱۷62۲ 0۷۵۲۵ ۴۲ 

6 ۲۱۵۷۵00۶ [1)2ع1د[ -4 

5- ۳۱۷۲010210 ۵۹۵0896 ۴ 


باشند. خروجی‌های مدل در هر زیر حوضه محاسبه شده برای هر واحد 
هیدرولوژیکی. به روش وزنی به‌دست می‌آیند (۲۸). پس از معرفی 
نقشه ۳1۷ با دقت ۲۰ متر که از اداره آب منطقه‌ای استان اصفهان 
تهیه شده بود» منطقه به ۱۸ زیرحوضه تقسیم شد. در ادامه شبکه 
آبراهه ترسیم و موقعیت دو ایستگاه هیدرومتری تنگ زردآلو واقع در 
نزدیکی خروجی حوزه آبخیز و پل بردکان در پایین دست سد و 
هم‌چنین موقعیت سد مخزنی قره آقاچ به مدل معرفی و مرز حوزه 
آبخیز تعریف شد. در مراحل بعدی» نقشه کاربری اراضی با استفاده از 
الگوریتم طبقه‌بندی غیرپارامتریک ماشین بردار پشتیبان (/53/۸ 4 در 
محیط نرم افزار ۷15.3[ تهیه شد. بر این اساس هفت نوع کاربری 
راضی برای حوزه آبخیز مهرگرد در نظر گرفته شد که شامل: کشت 
آبی (,[۸6/)» کشت دیم (۷۷۷/۲۲۲) باغ (016» سطوح آبی 
(۷۸۲» مرتع (0( اراضی سنگلاخی (406016 اراضی 
رشب[ 600 و خاک بانج له یه کین کتهانی سل 
معرفی شد. به دلیل عدم وجود نقشه خاک حوضه مورد نظر از نقشه 
خاک جهانی موجود در سایت مدل ٩۷۷۸۵۲‏ به شانی 
200 استفاده شد. هم‌چنین نقشه شیب حوضه 
بر اساس روش ۹05 به چهار طبقه تقسیم شد. در نهایت از تلفیق 
این سه لایه» نقشه واحدهای پاسخ هیدرولوژیکی (/1151) به دست 
آمدند. هم‌چنین حد آستانه برای کاربری اراضی ۱۰ درصدء خاک ۱۰ 
درصد و شیب ۵ درصد در نظر گرفته شد. بر این اساس ۵۴ واحد پاسخج 
هیدرولوژیکی برای کل حوضه به‌دست آمد. با استفاده از این حدود 
آستافر وی عاکرو کار راد فاطیی ‏ -خاگ و کاری 
اراضی در آن‌ها کمتر از حد آستانه تعریف شده است توسط یک روش 
میانگین گیری وزنی به خاک و کاربری اراضی که بیشترین مساحت را 
در آن حوضه دارند تبدیل می‌شوند. بعد از آن‌که از صحت داده‌های 
ورودی به مدل اطمینان حاصل شد نوبت معرفی اطلاعات اقلیمی به 
مدل است. در همین راستا با مراجعه به سازمان هواشناسی کل کشور 
اطلاعات ایستگاه‌های هواشناسی مورد استفاده مدل 518۸۲ برای 
بازه زمانی ۲۰۰۲ تا ۲۰۱۶ بر مبنای سال تأسیس و موجود بودن 
داده‌ها تهیه شد (جدول ۱). بر اين اساس بارش روزانه. حداقل و 
حداکثر دماء سرعت باده رطوبت نسبی و تابنش خورشیدی ایستگاه‌های 
موجود در منطقه برای بازه زمانی مذکور در قالب فرمت ۰16۲ وارد 
مدل شد. همچنین اطلاعات دو ایستگاه سینوپتیک شهرضا و بروجن 
به فایل ۷056۲ .۷/260 موجود در پایگاه داده مدل معرفی شد. هدف از 
انحام این مرحله استفاده از مولد داده اقلیمی ((۲۷۷07) برای 


بازسازی داده‌های مفقود می‌باشد. 


مصنط۱۷]۵ ۵۵۲۵۲ ۷ 5۱۵۵0۲۲ -6 
1 21067 ۱۷۷ -7 


اثر تغییرات اقلیمی بر رواناب و تبخیر و تعرق حوزه آبخیز مهرگرد ۲۲٩‏ 


جدول ۱- مشخصات ایستگاه‌های مورد استفاده در مدل 5۲۷۷۸۲ 
۱۵۵6 ۲۷۲۵۲ ع) جر 0ععته مصمتاهای عط) ۵ ی‌نامنه] 22 1۳6 -1 12016 


نام سازمان ارتفاع ایستگاه (متر) عرض جغرافیایی طول جغرافیایی نوع ایستگاه نام ایستگاه ردیف 
عصهه ممتاه 7 صهع ۵ رحعاعصص )طعتعط ممتاهاه تا ۱۰2 ت۱۱ ۹ عصهه ممناداه 10۳ 
مت 270 "31-25 510-47 سینوپتیک هیر 1 
00020( طافزه ۱۰۵ رصم 5 
دوانتانیی 1845 9-*31 "51-0 سینوپتیک ور 
/ ۱/۵2 0 ۵( 
وتان 2197 "31-7 51-87 4 ِ_ 3 
00020 2( 0 جهزر0ظ 
هواشناسی 1580 31-1 509-9 سینوپتیک لردگان 4 
00020۲( طافزه ۱۰۵ 0220[ 
و 2 2270 310-401 51-5 پارانتترتی اسفرجان 5 
20002( 6 1210[ ۱ 
وزارت نیرو _ 2130 31-0 510-0 بارانسنجی ری 6 
۵ ۴ ۷۲11507[ ماو 2[ تحطهو معصع [ 
وزارت نیرو _ 2230 "310-4 510-17 تبخیرسنجحی 0 7 
۵ 0 0۲ونط۷]1] 6۵ ۲۱۷۲۵۵0۲۵08 ۱۵۱۱۱ 
وزارت نیرو 2310 -*310 "510-447 تبخیرسنجی حنا 8 
۵ ۴ ۷۲1150۲[ مج ۲۱۷۲۵۵0۲۵0۲ ۲122 
٩0021 ۱/22۱7‏ ۵۷1 210 
وزارت رو 2270 "3216-29 510-87 هیدرومتری تنگ زردالو 10 
۵ 0 افنط] ۷ ۱۹۸ 0 20826 ]1 


به‌منظور براورد تبخیر و تعرق از روش هار گریوز-سامانی استفاده 
شد. هم‌چنین دبی ماهانه خروجی از دو خط لوله انتقال آب در حوضه 
و سد مخزنی قره آقاچ پس از انجام محاسبات اولیه به همراه 
خصوصیات فیزیکی سد به مدل معرفی شد. در نهایت مدل 5۷۷۸۲ 
برای اجرا آماده شد. 


واسنجی و اعتبارسنجی مدل 5۷۷۸1 

منظور از واسنجی مدلء بهینه‌سازی پارامترهای تأثیرگذار مدل 
است به‌طوری که خروجی‌های شبیه‌سازی شده بتواند روند و تغییرات 
داده‌های مشاهداتی را در بازه مورد نظر توجیه کند. اعتبارسنجی 
بررسی صحت مدل واسنجی شده در بازه زمانی دیگر می‌باشد 
به‌طوری که در آن متغیرهای شبیه‌سازی شده و مشاهداتی بدون 
تصحیح هیچ پارامتری مقایسه شوند. برنامه‌هایی به این منظور توسعه 
پیدا کرده‌اند که می‌تواند به‌طور اتوماتیک این عمل را در حوزه‌های 
آبخیز بسیار بزرگ انجام دهند. الگوریتم 511۳172 در بسته نرم‌افزاری 
٩۷۷۵۸ ۲-۳‏ یکی از اين برنامه‌ها می‌باشد که در منابع بسیاری از 
آن استفاده شده است. واسنجی و تعیین میزان عدم قطعیت. توسط 
معیارهایی به نام ۳-۳۵6۱0۲ و ۳۵010۲- بررسی می‌شود (۱). در این 
پژوهش با توجه به اینکه دو سال برای آموزش مدل در نظرگرفته شد 
بنابراین برای واسنحی مدل از بازه زمانی ۲۰۰۴ تا ۲۰۱۲ و برای 
اعتبارسنحی از بازه زمانی ۲۰۱۳ تا ۲۰۱۶ مربوط به ایستگاه 


هیدرومتری تنگ زردآلو واقع در نزدیکی خروجی حوزه آبخیز استفاده 
شد. در ابتدای اجرای مدل به منظور مشخص کردن پارامترهای موّثر 
بر رواناب منطقه» آنالیز حساسیت با استفاده از روش هر بار یک پارامتر 
(0۸) برای مدل انجام شد و پارامترهای حساس برای واسنجی 
مدل تعیین شد. در این پژوهش تعداد داده‌های مشاهداتی در خروجی 
به‌صورت ماهانه برای رواناب در مرحله واسنجی ۱۰۸ و در مرحله 
اعتبارسنحی ۴۸ نمونه مورد استفاده قرارگرفت. برای واسنحی مدل 
ابندا ۱۹ پارامتر که بیشترین تأثیر را در شبیه‌سازی دبی حوزه آبخیز 
مهرگرد داشتند انتخاب و بعد از انجام تحلیل حساسیت» ۱۰ پارامتر به 
عنوان تأثیرگذارترین پارامترها به مدل معرفی شدند. برای این کار از 
الگوریتم 51111-2 استفاده شد. 

الگوریتم 5171712 یک روش بهینه‌سازی معکوس است که درآن 
تمام منابع عدم قطعیت در دامنه معرفی شده برای هر پارامتر لحاظ 
می‌گردد. به‌نحوی که محدوده اولیه پارامترها بعد از هر بار تکرار 
به‌وسیله محدوده جدید به‌دست آمده توسط مدل, جایگزین شده و یک 
دامنه عدم قطعیت محدودتری از پارامترها به‌دست می‌آید. اين روند تا 
زمانی که دو شاخص کمی‌سازی عدم قطعیت به حد مطلوب خود 
برسند یعنی ۳-۳2610۲ به یک و ۳۵60۲-؟ به صفر میل کند ادامه 
پیدا می‌کند (۲۴). ولی اکثرا به دلیل خطاهای گوناگون در داده‌ها و 
عدم قطعیت مدل‌های مفهومی» هیچگاه اين مقادیر برای دو پارامتر 
ذکر شده به‌دست نمی‌آید. بدین منظور مدل مجددا اجرا شد تا بازه 
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مناسب برای هر پارامتر ب‌دست آمد. به منظور ارزیابی عملکرد مدل 
در مراحل واسنجی و اعتبارسنحی, از ضریب 15 ۵15 ۳۳۵6۲0۲ و 
1-6۲ استفاده شد. 


با انجام سعی و خطاء بهترین متغیرها برای کوچک مقیاس کردن 
متفیرهای دما و بارش حوزه آبخیز مهرگرد انتخاب و به مدل معرفی 
شد (جدول ۲). پس از بازتولید و آنالیز داده‌های تولید شده توسط مدل 
برای ارزیابی توانایی مدل از معیارهای آماری مانند 122 
5( و 5ل( استفاده شد (جدول ۳). 


شبیه‌سازی اقلیم آینده با استفاده از مدل 5۳5۲ 


صع) 0 صمااهازرنء تس 0۴ عصنایف‌مط سمل 10۲ ممی‌نلع صمناعع5۵ -2 12016 


جدول ۲- پیش بینی کننده‌های منتخب جهت ریزمقیاس‌نمایی بارش و دما 


جزئی اماره 0 مفیهوم پیش بینی کننده 
۳۵۳2 ۵ ملامناهماه ۳۲۵0۱0۱۷ 
۱۱ 
فشار متوسط یا 
01 001 روط نم برد 
16۷۵1 262 210 0۲6859۷۲۵ ۸۷۵۲۵۵۵ 
باد نصف النهاری 
0.068 0.02 ۱ 9 
۵ ۷611010021[ 
01 001 باد مداری ترا 103 ۱ 
50 ]مصصع‌نله آمازطاجه ۱۷۷10 
اگرایی تراز 02 850 
00۵4 03 با دی ۱ 
50 11۷01820000 احمحصصع تلظ 
ارتفاع ژتوپتانسیلی تراز 500102 
009 02 ۲صعصصع‌نله طعنمط ادتاهمام0۵ع0 
1802500 
04 00 باد مداری بر 102 ۳0 ۱ 
500 ]حمحصصع نله اهزاین ۱۷۷۱20 
ارتفاع ژتوپتانسیلی تراز 50010۵ 
02 015 ۲صعصمع‌ناه طعنعط امتامماه(مع0 
102500 
018 03 جهت باد ترا 0 ۱ 
2500 )جمحصصع نله ممتامعتل ۱۷۷100 
اگرایی تراز 02 
05 0006 ۵ 
0 11۷0۲2000 احمحصصع تلخ 
016 000 ی 
2850 ]حمحصصع‌تله مزا ۷۷20 
میانگین دمای تراز 2 متر 
02 021 2 اممصصصوناه معطهتهموصها وعمته۷ظ 
عنام 
03 000 توار متوسط فیح ریا 
۷۵ 162 21 0۲6550۲6 ۸۷۵۲226 
باد مداری تراز 002 
004 002 ۱۳ 
۱0 درمحصصع نله [ه)۶01ه ۱۷۷1۵۵ 
02 013 رتفاع ژئوتانسیلی ترا 
2 500 
04 05 باد نصف لهاری وت ۱ 
2صصم ناه 7/1۵0 همم ترع]۷( 
003 004 " رطوبت نسبی 
۷ ۱۵12101۷8 
میانگین دمای تراز ۲ متر 
01 05 


2 اصمصصصو‌تله ععتناموموها مممتو۷ظ 
عنام 


ته نوشت ری هو و رف 3 
کو 1 4 متغیر پیش‌بینی کننده متغیر 
۵ ۳۲۵/۱۵۱۷۲ ۲۰( 
۹30 
۱/۳30 2 000۳۹0 
۳۷ 76601-221 
28-0 11060008-6 بارش 
جمتاه]ز0 ]۳۲۲۵۵ 
8921 200002۵1 
830 166095002121 
۱۵۱۲ 11 _ 1060۳05 
8300 21.21 166۳09500 
۳50 0211 006005 حداقل دما 
حطاصتصز ]۷[ 
۳52۲ 26600528211 20 
۱ 11060008-6 
مصم 211 6۳0۵]) 1000 
۱/0 2 6۹10( 
۱۵۱۲ 11 _ 106005 
500 50216 
0 17060۹5052 بانط 
8 ی ات۷۲2[ 
۳ ی و 3 
ناه 1[ صصیاطوو260 
مصما 10606۳0۵2121 
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شکل ۲- مقایسه میانگین و انحراف معیار بارش حداقل و حداکثر دما مقادیر مشاهداتی و شبیه‌سازی 
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با توجه به نتایج موجود در جدول‌های ۲ و ۳ مشخص شد که 
داده‌های درجه حرارت همبستگی بهتری با داده‌های مشاهداتی در 
مقایسه با داده‌های بارندگی داشته‌اند. به عبارتی تغییرپذیری دما نسبت 
به بارندگی کمتر است و دما یک متغیر پیوسته است که کمتر تحت 
تآثیر ناهنجاری‌های زمانی قرار می‌گیرد. درحالی که بارش تحت تأثیر 
عوامل مختلفی در حوضه قرار می‌گیرد. از سوی دیگر, بارندگی یک 
متغیر گسسته است. نتایج بسیاری از مدل‌های گردش عمومی جو و 


مدل‌های منطقه‌ای نشان می‌دهد که مدل‌سازی بارش معمولا همراه 
با خطا صورت می‌گیرد. هم‌چنین مدل ٩19۷‏ یک مدل 
رگرسیونی است. این مدل‌ها اکثرا با ایجاد روابط رگرسیونی بین 
متغیرهای بزرگ مقیاس و ریزمقیاس تعریف می‌شوند. برای 
شبیه‌سازی بارش» ریزمقیاس‌نمایی وقایع حداکثر در آینده توسط 
مدل‌های رگرسیونی تا حدی مشکل است (۳۵). لازم به ذکر است که 
حوزه آبخیز مهرگرده کوهستانی است بنابراین بارش را به شدت تحت 


۳ شریه آب و خاک. جلد ۲۴ شماره ۱. فروردین - اردیبهشت ۱۳۹۹ 


تأثیر قرار می‌دهد و احتمال بروز خطا را در اين متغیر افزایش می‌دهد. 
به بیان دیگر مدل ٩۳51‏ نسبت به شبیه‌سازی بارش از حساسیت 
بالایی برخوردار بوده لذا هر عاملی که اين متغیر را تحت تاثیر قرار 
دهد امکان بروز خطا توسط این مدل را تشدید می‌کند. به همین علت 
شاخص آماری نش ساتکلیف در ارتباط با بارش نسبت به حداقل و 
حداکثر دما پایین‌تر برآورد شده است. که با توجه به ماهیت مدل و 
شرایط توپوگرافی حوزه آبخیز این حد از اختلاف قابل قبول است. در 
این راستا می‌توان به مطالعاتی چون عساکره و اکبرزاده (ه) در 
ایستگاه سینوپتیک تبریز» نادری و همکاران (۲۷) در حوضه سیمره 
اشاره کرد. بر این اساس شاخص‌های آماری به‌دست آمده برای حوزه 
آبخیز مهرگرد در سطح قابل قبولی قرار دارند. در شکل ۲ به خوبی 
واضح است که میزان بارش شبیه‌سازی شده برای آینده نسبت به 
میزان بارش مشاهداتی به‌طور قابل توجهی کاهش می‌بابد. هم‌چنین 
بیشترین اختلاف در ارتباط با انحراف معیار بارش مربوط به ماه‌های 
ژانویه و دسامبر به ترتیب به میزان ۴/۱۴ و ۴/۳۷ واحد کمتر از مقدار 
انحراف معیار مقادیر مشاهداتی است. در ارتباط با دماء که در برخی از 
ماه‌هاء خصوصا در ماه‌های ژانویه و دسامبر درمورد دمای حداقل و در 
ماه‌های ژانوبه و نوامبر در مورد دمای حداکثر اختلاف نسبتا کمی 


بین انحراف معیار مقادیر مدل‌سازی شده و مشاهداتی وجود دارد با 
این وجود مدل‌سازی دماهای حداقل و حداکثر به میزان بسیار خوبی با 
داده‌های مشاهداتی مطابقت دارد. نتایج ریزمقیاس‌نمایی مدل گردش 
عمومی جو در این حوزه آبخیز نشان‌دهنده کاهش ۵۳/۴۸ درصدی 
بارندگی و افزايش حداقل و حداکثر دما به ترتیب به میزان ۰/۸۴ و 
۹ درجه سانتی گراد در دوره شبیه‌سازی شده ۲۰۳۰ -۲۰۱۷ نسبت 


به گذشته است. 


شبیه‌سازی شرایط هیدرولوژی منطقه با استفاده از مدل 
۸ 

با توجه به این‌که ایستگاه هیدرومتری تنگ زردآلو به لحاظ 
کیفیت و طول دوره آماری در حد مطلوبی قرار دارد بنابراین از 
داده‌های این ایستگاه برای واسنجی و اعتبارسنجی رواناب ماهانه 
استفاده شد. بر اساس نتایج حاصل از آنالیز حساسیت مدلء ۱۰ پارامتر 
به‌عنوان حساس‌ترین پارامترها شناسایی شدند. که نتایج آن در جدول 
۴ ارائه شده است. هم‌چنین نمودار همبستگی رواناب مشاهداتی و 
شبیه‌سازی شده مربوط به دوره واسنجی و اعتبارسنجی ایستگاه تنگ 
زردالو فر شکل ۲ ارائه شنده است, 
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شکل ۳- نتایج واسنجی و اعتبارسنجی رواناب ماهانه ایستگاه تنگ زردآلو با استفاده از الگوریتم 5111-2 
صطاتمواه 2-]۱۳ 6ظ۳) عصنعه ممتاهای وماق2270 عمج ۱ صد امد رواجم ما ۵ المع ممتاه0ناه۲ مه ممتاهرطانامن -3 مساهز۲] 


با توجه به نتایج موجود در شکل ۳ شاخص ۳-۳۵0۲0۲ -؟1 
۲ 1 و ۸5 برای دوره واسنجی ۰/۵۲ ۰/۲۴ ۰/۷۲ و ۰/۶۹ و 
برای دوره اعتبارسنجی ۰/۴۵ ۰/۲۹ ۰/۷۱ و ۰/۵۸ به‌دست آمد. 
بنابراین درصد قابل قبولی از داده‌های مشاهداتی در باند عدم قطعیت 
۵ درصد قرارگرفتند. هم‌چنین نتایج نشان‌دهنده این است که مدل 
به‌صورت رضایت بخشی رواناب ماهانه حوزه آبخیز مهرگرد را 
شبیه‌سازی کرده است. با مقایسه هیدروگراف شبیه‌سازی شده و 
مشاهداتی به خوبی واضح است که در بعضی از ماه‌های سال (می 
۴ و مارس ۲۰۰۶) برای دوره واسنجی و (می ۲۰۱۳ 
مارس ۰۲۰۱۴ آوریل ۲۰۱۵ و مارس۲۰۱۶) برای دوره اعتبارسنجی در 
قاط پیک, تا حدودی همبستگی کمتری دارند. در اين ارتباط می‌توان 
به مطالعات تولسون و شوماکر (۳۳) و چو و شیرمحمدی (۷) اشاره 
کرد. دلایل این مسئله عبارت است از: کوهستانی بودن و وجود سنگ 
آهک. آهک مارنی و آهک دولومیتی مربوط به کرتاسه زیرین که از 
ضخامت قابل توجهی برخوردار بوده و بخش عمده نواحی مرتفع را 
تشکیل می‌دهند. بخش‌هایی از جریان که مدل با اختلاف کمی نسبت 
به هیدروگراف مشاهداتی شبیه‌سازی کرده است. در اواخر زمستان و 
اوایل بهار به خصوص در ماه‌های (مارس» آوریل و می) اتفاتق افتاده 
است. از عواملی که سبب افت هیدروگراف مربوط به دبی شبیه‌سازی 


شده است» برداشت آب با هدف تأمین آب شرب مناطق مجاور 
همچون شهرضا و سمیرم با استفاده از خطوط انتقال آب و احداث سد 
مخزنی قره‌آقاچم جهت تأمین مصارف مختلف ساکنین حوزه آبخیز 
می‌باشد. 


بررسی آثر تغییر اقلیم بر رواناب و تبخیر و تعرق 

با توجه به شکل ۴ نتایج به‌دست آمده از شبیه‌سازی اقلیم برای 
بازه زمانی ۲۰۱۷ تا ۲۰۳۰ توسط مدل 55۷ حاکی از کاهش 
بارندگی و افزایش دما در تمام ماه‌های سال در مقایسه با دوره پایه 
می‌باشد. بر این اساس شبیه‌سازی رواناب حوزه آبخیز در دوره آتی و 
سقایبة مقاذیر آن نا فوزه بایه ببانگر کاهتن شورس رولاب مالاند 
درازمدت در دوره آنی به میزان ۲۳/۸۲ درصد می‌باشد. هم‌چنین پیک 
دلیل شرایط دمایی گذشته بیشترین میزان رواناب در اوایل فصل بهار 
همزمان با ذوب برف زمستانه و بارش‌های بهاری اتفاق افتاده است. 
درحالی‌که در آینده نزدیک هرچند بارش‌ها از نظر کمی نسبت به 
گذشته کاهش می‌یابد ولی به دلیل تغییر رژیم بارش‌ها» نقطه اوج 
رواناب از اوایل بهار به سمت اواخر زمستان منتقل می‌شود. این نتایج 


۴ نشربه آب و خاک. جلد ۰۳۴ شماره . فروردین - اردیبهشت ۱۳۹۹ 


با مطالعات منصوری و همکاران (۲۱) در حوضه زرینه رود اکبری و 
همکاران (۳) در حوضه رودخانه ارس و امیل سورو و همکاران )٩(‏ در 
حوضه 52002۳12 مطابقت دارد. با توجه به اينکه تغییرات دمایی 
عامل اصلی تغییرات اقلیمی شناخته شده است. معادله‌هایی که بر پایه 
پارامترهای دمایی باشند. با توجه به پدیده گرمایش جهانی و تغییرات 
اقلیمی قادر به پیش‌بینی میزان تبخیر و تعرق در دوره‌های آتی خواهد 
بود لذا از روش هارگریوز-سامانی که روش دمایی است به دلیل 
نیازمند بودن به داده‌های ورودی کمتر و پاسخگو بودن به مسئله 
پدیده گرمایش جهانی و محدود بودن داده‌های ریزمقیاس شده در 
حوضه برای بررسی اثرات تغییرات اقلیمی بر تبخیر و تعرق استفاده 
شد (۱۱). بر این اساس نتایج نشان‌دهنده کاهش تبخیر و تعرق واقعی 
آینده ۲۰۱۷ تا ۲۰۳۰ در مقایسه با دوره پایه ۲۰۰۴ تا ۲۰۱۶ است. 
به‌طورکلی نتایج نشان‌دهنده کاهش متوسط سالیانه تبخیر و تعرق 
واقعی آینده در مقایسه با دوره پایه به میزان ۲۶/۰۳ درصد می‌باشد. 
این مسئله ناشی از کاهش بارش در آینده و عدم وجود آب کافی به 
جهت انجام مکانیسم تبخیر و تعرق واقعی می‌باشد. بنابراین با توجه 
به بحث افزايش دما و کاهش بارندگی برای دوره آینده و به دنبال آن 
کاهش رواناب سطحی و زیرسطحی این مسئله خود کاهش رطوبت 
خاک و در نتیجه کاهش مکانیسم تبخیر و تعرق واقعی را موجب 
خواهد شد. با توجه به شرایط اقلیمی حوزه آبخیز مهرگرد بیشتر 
بارش‌ها در فصل زمستان به خصوص دی و بهمن ماه اتفاق می‌افتد و 
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در فصل تابستان از میزان 
نحوی که کمترین میزان بارش مربوط به ماه‌های تیر و مرداد است. 
این حوزه آبخیز بر اساس تقسیم‌بندی کوپن جزء مناطق معتدل سرد با 
تابستان‌های گرم و خشک و در ارتفاعات جنوبی جزء مناطق معتدل 
سرد با تابستان‌های خنک و خشک می‌باشد. همچنین رشد پوشش 


بارش‌ها به شدت کاسته خواهد شد به 


گیاهی در اين حوزه آبخیز از اواخر اسفند ماه شروع و در اواسط 
اردیبهشت به حداکثر رشد رویشی خود می‌رسد (۰۲ ۱۸ و ۲۰ بر این 
اساس بیشترین میزان تبخیر و تعرق واقعی در این حوزه آبخیز مربوط 
به اواخر زمستان تا اواسط به‌خصوص ماه‌های فروردین و 
اردیبهشت می‌باشد. علت این مسئله رسیدن پوشش گیاهی به حداکثر 
رشد و وجود آب کافی است. در ارتباط با بحث تبخیر و تعرق پتانسیل 
نیز باید گفت با توجه به ارتباط تنگاتنگی که میان تبخیر و تعرق 
پتانسیل و دما وجود دارد بنابراین تبخیر و تعرق پتانسیل سالانه نسبت 
به دوره پایه به میزان ۱۰/۲۰ درصد افزايش می‌یابد. بر این اساس 
میزان تیحیر و تفرق پتاسیل در حوزه یز مهرگرد در فسل تایبعان 
با افزايش دما به حداکثر خود خواهد رسید. . نتایج به‌دست آمده با 
مطالعات سیفی و همکاران (۲۲) در ایستگاه کرمان» کوهی و ثنایی 
نژاد (۱۹) در حوضه رودخانه ارومیه. گودرزی و همکاران (۱۲) در 
حوضه کرخه و لیانگ و همکاران (۲۰) در حوضه رودخانه 12067 در 


این باره هم‌سو بود. 
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شکل >- تغییرات رواناب » تبخیر و تعرق واقعی و پتانسیل تحت تأًثیر تغییر اقلیم (۲۰۱۷-۲۰۳۰:۲۰۰-۲۰۱) 
(2017-2030 ,2004-2016) 2086 عاقصتل ۲ 216060 اهتاصهامم مه مهتم تمه موی ماه معصمت م۲2 10۵1۴ -4 موز 


اثر تغییرات اقلیمی بر رواناب و تبخیر و تعرق حوزه آبخیز مهرگرد ۲۳۵ 


نتیجه گیری به بارش اثرگذاری عوامل مختلف بر بارش هم‌چون توپوگرافی 


۰ ِ 2 توان گفت که مقدار به‌دست آمده برای ضریب نش -ساتکلیف قا 
ده تس می‌توان رب برای ضریب نش یف قابل 


اصلی بشر در دوره‌های آتی به حساب می‌آید جهت مدیریت هرچه 
بهتره ضرورت دارد وقوع تغییرات اقلیمی و میزان تآثیرگذاری آن بر 
منابع آب ارزیابی گردد. در این پژوهش سعی شد تا با استفاده از مدل 
اقلیمی 020۳5۸2 و مدل هیدرولوژیکی ۸۲2012 5۱۷ اثر تغییر 
متغیرهای اقلیمی در آینده بر بیلان آب حوزه آبخیز مهرگرد 
به‌خصوص رواناب و تبخیر و تعرق (واقعی و پتانسیل) مورد بررسی 
قرار گیرد. 

براساس خروجی‌های مدل 50584 مشخص شد که بارش در 
مقایسه با دوره مشاهداتی به میزان ۵۲/۴۸ درصد کاهش و دمای 
حداقل و حداکثر با اختلاف کمی نسبت به دوره مشاهداتی به ترتیب 
به میزان ۰/۸۴ و ۳/۹۹ درجه سانتی‌گراد افزايش می‌یابد. نتایج حاصل 
اضیا ها آمارف اک ار ان اسنت که موه ی های اقلیمی 
حوضه را به خوبی شبیه‌سازی کرده است. بر این اساس ضریب 15 
2 و 38 برای بارش به ترتیب معادل ۰/۹۲ ۵/۸۱ و ۰/۳۹ و 
برای حداقل و حداکثر دما به ترتیب معادل ۰/۹۹ ۰۰/۱۶ ۰/۹ و 
۹ ۰/۲۱ ۰/۹۹ به‌دست آمد. اگرچه ضریب نش-ساتکلیف در 
ارتباط با متغیر بارش در مقایسه با حداقل و حداکثر دماء پایین براورد 
شده است اما با توجه به ماهیت مدل ]1291۷ حساسیت مدل نسبت 


قبول است. به بیان دیگر مدل 505 در شبی‌سازی بارش از دقمت 
کمتری برخوردار است. در بخش هیدرولوژیکی معیارهای آماری "13 
ول ۳-۳2۵۱0۲ و 3-۳۲260۲0۲ برای دوره واسنجی به ترتیب ۰۰/۷۳ 
۵ ۰/۵۲و ۰/۲۴ و برای دوره اعتبارسنجی به ترئیب ۰/۵۸۰/۷۱ 
۵ ۰/۲۹ به دست آمد. لازم به ذکر است که مدل 5۷۷۸۵۸۲ 
بارش را با استفاده از متوسط دمای روزانه و به‌صورت باران يا برف 


تقسیم‌بندی می‌کند. با توجه به ساز و کار مدل در ارتباط با شبیه‌سازی 
رواناب ماهانه ماه‌های پس از فصل سرما که در طی آن‌ها با گرم 
شدن هواء برف باریده شده در فصل سرماء ذوب می‌شوده می‌توان 
نتیجه گرفت که مدل قادر به شبیه‌سازی خوب ذوب برف در برخی از 
تامهای میا هه اش مود خربای ای یدبا رش قیده فطاق: 
کمتری با مقادیر مشاهداتی دارد. در ادامه با هدف بررسی اثر تغییر 
اقلیم بر رواناب و تبخیر و تعرق واقعی و پتانسیل حوضه اقلیم شبیه 
سازی شده برای آینده به مدل ٩۷۷۸۲‏ معرفی شد. پس از تحلیل 
نتایج مشخص شد که متوسط رواناب سالانه به میزان ۲۳/۸۲ درصد 
روند رو به کاهشی داشته است. هم‌چنین تبخیر و تعرق واقعی به 


میزان ۲۶/۰۳ درصد کاهش و تبخیر و تعرق پتانسیل به میزان ۱۰/۲۰ 
درصد افزايیش یافته است. 
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79-۰ :326 ۲۲۵۲۵۱۵۵۷ ۵۶ (مصیامل .ماطصتامی طفتاترظ رصتعهه صمعصقان مطا ۶ه و0عطوعته/« حاباو صا معصقطه متفه 

۲ :5 ما۷0 .صتفوظ صناتفک کجعرن) قفقعته ناه ۵۶ معصملده قعم‌تناموع۲ ۳/2 0۴ ]1586001 .2017 ۲0۷96۲۰ ۵۶ «افتصتا -24 
-105-> :0۵06 .(2324) هعتظ لیگ 0عموعطع]۱۷ مطا ۵۶ عمصفلهه عمج ۷۷ :24 جرخ ]۵۵0۲ اصمحصووهووظ وععتا0وم۴ 
.39 

0 .۸ ۷00]۵62010۵4- ۱۷۲۵212۷ و.ل ۵62-۳۱۵۱۵۵48062 ریگ ۱۷2۳062-6۲۵2 ری هزتوم۱ ۲۱۵1۱201۵ ,.ظ مهد مصا۱۷]0 -25 
۰ ۵271و آهعنانهمتبوه مه صا وامدمرحصا مفعصقطه مقصتا مه مصتام0مصظ لهمنوم۳۱۲0۲۵۱ .2015 هر طع۷-عتاوه 
29-۰ :175 مممصه۱۷]2 ماه ۷۷ .مخ .۱۷۲۵۵0۵ رصتفهه ععبنک ممتالممنان :تیاه وققت) 

۰ )هنن تمه م2 ۳۱۷۵۲۵1۵۵10 .2012 .ظ 00۱۷۵2 20 ر.ل ۸010 ریک صتاصه-عمامداهنا ر.رظ ممولز ۷۷ .نا تففت:1۷0 -26 
.1241-7۰ :(55)4 ۲ظهمف .عصه 1 بصمت)اع۷211 20 محمتاه‌طلامه رمونا 

متقصره5 طاً دتمامصصه‌تدم عاقصصتان مم معصقطاه ماقصصتن 0۴ ۴۲۴۵۵۱5 .2018 ,۱۷۲ مصجامنا تصصهل‌هطن 24 .۷ ۵۵02۲21 رو ت20 -27 
ماه طوتامدط! طازه صدتوتعظ 18) .69-76 :(11)039 عصارمعصاعصظ ی ممصمزهد اصمصم‌مهم۱۷۸ 0عطوع]2 1۳2۵-۷۷ .صتفح‌ظ 

,20101) رصاعم وم امممصرووعووه همه 20 [زمو .2005 بکز.ژ فصصهیللز۷۷ 20 ریک.ل ننک ریت۲ ۳010 ولو واه -286 
۰ 494 .(11201121 5" 12961۲) 

0۵1 :0۵01 )همصصوومووه ماه 20 لزمو ,2011 .ل فصصهتالز ۱۷۷ 220 و.ل تتصنک ر.ل 10مصه 20 ریق طعفالع -29 
۲۰ ,)تافص موم م26 ۷۷ عق1۵ ,2009 مملوع۷ بطملاه)مممصتهم 

قیامصمعتصا "قعصصمل ,2019 به لمامعقصصطع 1 جرج یه تصفامع7 92 تلعهاه ۳ تامطومة‌مطل؟ ربا تصعوظ بب تطممطح9 -30 
0 224 ۲۳60۲۵621 .صوع1 ل2اجمن صا جات مدع صا فععصقطه متقصصاه ما صمتاه)مم0ه تتقطا مه 1۳20۷716026 
1377-2-4 :137 7«و0۱0)مصن) 

۰ 7/۵۱۵ ۱۵۷۵۱ متفجه بم۶ ام0مصه 5۱۷۷۸۵۸۲ 0ع)2ععاص1 ورن ۶ه ممتامه‌تانم۸ .2009 .ظ مامصصعبوطه مه ریکانکا صفهتطاهه -31 
100-۰ :37 ۵۵86۲۷2008 5011 گه اقصتامژ ق1801 

کلو مه 12010۲ 20 طملووع۲وع1 عدعصناً 0۶ ۵061 0عصنطاصم 2 عممم12۵۷۵۱ .2011 .۳ تطوتک مه ,9 ما۱۷ ره طلمه5 -32 
11 هه ععند۷۷ ۵۶ فتاه[ .رصمتاهاگ صقصصعفک :واگ معمت) ممتاهتامومماممد۵ ممصم۲۵۶6 ممتامنام7ه 10۲ (هعم-1۲) 
رامجتدطاه مفتاعصط طازه۳ صهتوع۴ 12) .1186-1196 5016806 

۵61 :5۱۷۷۸۰۲2000 وتا صتعوه ع111بصمرص فط ۵۶ عصتاممصص 0مصفتماع ۷۷ .2004 هی عملقمصهمطو مه رف.ظ صمقام1 -33 
فطا ما ۲ممفصه فمدام‌ومطن مج تممصصتلعی سمل ۵ ممتا‌تلمتج مط مر ممتلهل‌ناه۷ مه ممتاه‌طاتلهه رتصعصصمماهبع0 
۰ 159 .عصازممصلمص مامممممبصو مه زره ۵۶ اممطمو ,۲۵9۵۲۷01۲ ۷1۱1 فحمصصوم 

ماه عم فتورلممه بامتدهعص مه صمتاهه‌تاصعل1 تعامصح‌تعظ .2015 اه قمع ۷ مصه رنکا,۷ فطل ر.ظ لهتصنا -34 
.3744-6۰ :(29)017 ۲۶۵۵۵9969 ۳۱۷0۲0۱021021 رهتلصا صعاعحظ ۶ صتفدها بع ۲۱۷ ۵ صا ماهتا 
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آمصمنعع۲ ۵۶ ]همصرووعوعه مطا عم اما ٩۱۵۵0۲۲‏ 06019101 ۸-]9۷نا .2002 ۴.۷۰ 0۷تتدظ مه .2۰۷۷ 120۷908 وینآنک ۷۷110۷ -35 
٩010۷۵۲۵ 17: 147-۰.‏ ک مصنام0مجظ لهاصمصصممت هو رواعدمرحصا معصقطه متفه 

متقصصتاه لقجمنوع: ۵۶ )ممصصووعووه فط 10۲ ۵01 ۹۱۵۵۵۵۲۲ 06019101۲ ۸ -4.2 9125۷۲ .2007 ۷۷۰۵۰ 12۱۷908 20 رانک ۷۷119 -36 
۰ 94 ,۷۷۵۱69 0ج 0صداعدظ ۵۶ رمجه‌عخ )صمصصدم1ص 4.2 ممتفتع۷ آمناصفجط ]9۷اه رقامهمصا معصقطه 

2 :همع ۴ لصا دسمل میاه مهم معصقطه متفصتا ۵۶ قاعجمرحص1 .2016 .1 0۵ 220 وت صقن ریم ۷۵۵ .۷ عطه2 -37 
۰ بر تهع دمک مرعمطامومصصخ نصا کههم2۵ 10 بقصتطی رصتفحظ هنک تک صا رتاو عقهم 

۷۵۲ ۰ ۷۷ ظا طا معصقطه متقصصتاه م6 ۴مصند ۵۶ مقممجوم ,2015 .۲ وصج۷ 20 ریک عناممجح‌طاهاض رل مقطم و2 تاو ,رن 200 -38 
37-۰ ل2ظا10 وعمممزهی لهعنع0۲010 ۲۱۲ .مصنطن رصتفهط 


نشریه آب و خاک (علوم و صنابع کشاورزی) انمک 0جد ۱۷2/۵ ۵۲ احسمز 
جلد ۰۳۴ شماره ۱. فروردین - اردیبهشت ۱۳۹۹ ص. ۲۲۵-۲۳۹ ۳۹| 225-9 0۰ ,2020 1۷۲2100-4011 با ۷0۰ و34 ۲۷۵0۲۰ 


۵ ۱2۱26 و۳۵۲۵۱۸۵۵)۵ ع)مصصنان م60 0۴ ۲۳۲۵۵ 0۴ وزورادحخض 
۵ ۱۷۲۵۱۲۵۵ 0 ۵۲۲۹۵۵۵۸ ۲۱۱۲۷۵۵0۵۵۹۵۱۵ 


تتصنطفرطظ 1 ۲2۱۵0 هل ,7 
02-09-09 :160661۷60 
12-01-2020 :660۲60 


فص 0۷۵۲ وملهعو ل2طاماع مه لمجمنوع طاهطا ۵ عصتع‌صمل معط مق متقصصتاه مطا توص ]فوم معط ما تممتامعت0 1۸0۳0 

۵ ۱06 ۲6 و1 معصقطن متقحصتات .عتنظ عقعص فط صاً میمتاومی آذه وععصفطا میاه قط 0ع۲ععهرم 0فل2 و1 16 .انوم 
معط 0 ۵۱2۲60 لمع ود وعوهع فعفط ۵۶ ممتاحتامععوهه فط1. ,عتعطام‌ومصاه فط صا فممتعوتصه ففع وفبامطصهعو 
۰6 6ظ) رتصنامحصح رمصصا صا فععصفهه طفامنط) مآمین لمعزوما0»0 بط مطا مامعله معصفطل میفصصنلت .عفهععمص متتلج همم 


۵05 ۷0۲0101621 0۶ ۵و۱ عظ) رسه100 .عاه مناد سوام 5۵ تمهت 20 ماه ممتهاممه۱ه رصمتاهاتم0ع ۵0 
۵۶۲ م20 ارم 1 .0عطوعمنه/< معط ما ملمه لهعتوم070۱ظ معط فصتامعله و1200 مط) 2۷۵ظ 60 06۷610060 صععه 2۷۵ 
۶ وا0عه فطا عصتووعوعه ۵۶ 0مطامصظ حرمحصجصرمی ف1. .قمع معط ۵۶ مامجصقره صه و1 (6۷۷۸۲) ۰۲001 )صعحطوو6ووخ 


2101 0۳ (00۷۲5) 0615مظ ممتاهانمیتن لحتمصعع از عصمله واعل‌مصه لمعلو1:0۱۵بوظ عصتوده ود سم وم ععصقطه متقصصئاه 
0 ۲مصنا جم فعصقطه مافصت 01 امعه عمط ماهعتاده1۳۷ ۵) و1 رتاو فنط) ۵۶ مو۵ماض مظ1۳ .(۴0۵) واع0مصه تمعطهع۷ 
عط لمح عم منوما00ظ 5۷۷۸۲ فطا عصلونه 0عطوماه ۷۷ موه ۵۶ (لقتاصعامم مه لهعع) متاه ام مممام هه 
1۰ 011502110 ۱۷۲2وصمن) 
0 ۲ عمط ما ومامصصهتوم متقصتاه آمممنوم: ۵۶ مله: معصجطه مطا هام2۵ ۳۵۶ :۷۲۵0۵09 200 واهزه۱۷2)6 
فصح تمه که اج لتق فطا رقتعامصومتهم لممنع۸0۲0۱۵ظ جوم ومعصقط ففقطا ۵۶ )0عه فطا 20 (2017-2030) 
0 0296 6ظ) 10۴ جمتاهاه م2272 معصه 1 مه ۵۶ ممتاهاتوهع۳ن 0صه صمتنهاه ممزنت0ظ مطا ۵۶ مجومصهه تال ت۳۵ 
عط ۵06۲ 9125۷۲5۰2 ۶ه 2۵061 مطا مصتوتا رراعصتلتهعع ,مه 9۷2صدن فطا ۵۶ قاتامرصا قح عون عت۵ (1984-2005) 
عم 5۱۷۲نا مطا ۶ بمصمتم‌تلله معط ملقلد7ع 10 بعمتاهعوه مصتاهعصومل ما منوج ۳2۵ ۴۵۲8.5 ۵۴ متقطععو 
عتقلنصته ما 1260 مهب« محر 2012 6۷۷۸۲ ما ,معاد دعه عطا ما .۱۱6 0صح ,رظ8] ,۳۳ متماتت لمه‌تامتاهاه ۵0و مه 
0 ۷۷2۶ ۲6۵10۲ عطا روعع‌اممط 20 ۵۶ مملوه0۳6 ۵ ۳/165 مقصه 1۱۱۳۷۲ مطا عصمن0هعصا ماخ .ععمتاتجمن منع۶۱0۲010ظ عطا 


6 1۲0۱0۵102 ۵۶ واتصت 534 ر1006٩‏ م2 باه ,دصق موب 0صق1 گ۶ه ممتامصتطاجم فط) مصمعر .عصلفعه-طاناه 18 1۳10 


6 6۲۵60۱۲۵ ۵ع] متاحل رم13 لصج مصیاحصتصتمه رطمتامازمتمع۳ مطلل تنامض 0242 مامتان مصفط 1 .0مصتهاهاه مععبه رتا۳۴ظ) 
6 20 صصحل معط ۵۶ ععهعوع۲ عط م60 هبات ,90061 عه م6 1000006۵60 ۷۸۵۲۵ ۲20120109 تقامو 20 رلععوو 0ه1 ,اتتصصباط 
1۰ 9۱6 میهد 006۵0 مهد ره همه م۱۷۵۲ ممههطمونل ره هه اهملووطام بهعته مطا ما دعصنا تعاعصهتا م۷۵ 07۷0 
۵ 0۴0668 0120108 له فظ1. ,حصطتنمماه ناو رها مصمل. مهس اعمصصر قط گم ممتامع‌تله۷ مه صمتاهءطانلم مط1" 
علونله7ه 10 0۲۳06۲ 18 .2016 1۵ 2013 ص۱0 ۷/28 ووععمتم ممتاهعناه۷ عظ) ملنط7 2012 6۵ 2004 ۵67100 عظ 10۶ 0ع)عنلههم 
01 معا رعدمم‌تم فنص ۲۵۲ .0ممنا متعه جماع۳-۲۵ 4صه ۲مام۳-۳۵ظ ۲۳ بواظ ملصعت‌ت1له0ه ,لمح معط آن ععصفصمهج عط 
۵۵ ۵061 ۳۷۲0۱021621 عطا عصتاهتطانلهه ماخ .عامصصهتهج طمقه ۲۵۲ معصهد ماقتعو0؟وم2 معط صتقاهام مه 9)27160ع1 ۷729 
عم لعاهعناوه 1۳۷ ۷۷۵۵ ععصقهه ماقصئله ۵۶ ]6۴۱60 مطا رتالمصز۲ .۵061 5۱۷۷۸۲ مطا مه متفصصتاه 0عهانبصته مط) 1۳0۲000660 
۰ ۱۷۱۵۲۵ ۵۶ (افتاممامن 0ج لهم) ممتاهامومهت0م ۵2۵ 20 ۵0۴صنا۲ 

2 18016206 ممتععر عتطا صا 0۵1مصر متقصناه مطا ۵۶ عصنلههعه0۵ فا ۵۶ فالناعع: ۲ خصمزقع کر م0صد والنعع1 
0 ماج 0۲ ومعتهتهم‌جصع) امک تمطظر لمح حصیاصتصتطظ . وففعتمص1 4صه صمتامالمم0:6 0 53.480 0۴ 060۲6286 
,۵061 5۱۷۷۸۲ عطا 0۶ ولورلمصه بدتاتوجعی معط ۵۶ فاناوع مط) وم 860ظ .۲۵5۵600۷۵ ,3.996 هه 0.84 ,(2017-2030) 
همان دمتادتاهاه مط) بممتامعه لهعزع1۲0۱0 ۲ معط ص .قتماممصه‌نهم 101۷۵فصعو ]دمص فطل و2 1060001160 ۷۷۵۲۵ وععامصصه‌نهم 10 
۷617 ,0.24 200 0.52 ,0.69 ,0.73 06700 صمتاحطنلی عطا ۲0۲ 0عصنهاداه ۵۵ ۲-۲۵0۲ 8صه عماع۳۲2-ظ روا ,۲2 آن 
بانط مطا ما ممتاهانصته امد مدوم بواه1۵50۵6۵0۷ ,0.29 0ص 0.45 ,0.98 م0۰71 ,06100 ۷۵11021108 6ظ) 10۶ 20 
٩00۷۷64 2‏ ۵100 عفهه عطا طز ممبله۷ فلا ۵۶ ممفتتهمرجم 20 معصفطه مفلقصتا ۵۶ معصمتاکصا عص مصنا ۵۲100 
۲ ممتاهانکصمتممجه تمتتامه نت۲60 مفله ]7۷۱1 ععصقطه متقصصات ,1۶مص ۵۶ ممدعهه تمیاصصه صه طاً 23.826 0۴ 06016286 
۰ 0 1۲۵)10۴مفصمتاه مج ۵ آمتاجعاه0 عفععتعص1 همه 26.030 


0 ,1665011665 لهعنهآظ گه بالنامو۲ رعطلهممتمط۴ )ممصمم ۷۱۸۵۵ 20۵۲۹۳60 ۷۷ صا ر0ووع۶۵1ظ جرج 60)ه20زن0 ,۷۲,96۵ -2 4صه 1 
۷۰ ,۲11۷۵۲۹16۷ 

(6۵۷۵270.260.11) 2160181 :اتمحطظ ماخ مصتمممومتمن -۶) 

«صجصصصمن) ماع ۷۷ آمطمتوم۴ صقطفکو1 رمع‌تلن ولو رلممض متصجمجمعظ قصه عصتصصهاظ عط) ۵۶ ۲۱۵2۸0 ررع۱0مطم۲مه۵ع۵) ط ,2.ظ -3 
9 :]10۵ 
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۲ و1 ممتاهاز۵ع0 عطا عقطا 0مصتجماعل وه ]1 رلفل‌مصط ]ناو فطا ۵۶ فالناوع: فطل رم 800ظ مصمتعسمومن) 
1۱06286 ومناه ومع متام هم صرح تنم مط) لمع رلمتوس مناج واه مه طز جموتتهم‌صمن صا ۵0منا۲۵0 
فقط ۵061 9۱۷۲5 مطا رهتتعاتت لهعتافتاهاه ما عصتلت0ععخ .مه ممتاح وان فط م6 ۵0عوم‌جمع عمصمممن طوتله 2 
۲ ۱۱5 مه ,۳۵/۹۴ ,2 ۶و معله معط رواع‌صنلتمعع۸ ,منهج عتععبط عطا عم متعتعطصفتمم فطا عصتاماصته مد 666064ع90و 
صتتطر لت حطظر. صرق مصیاطصتطتص مطا تم هه ۲6506۵ ,0.39 20 5.81 ,0.92 0 لهنا0ع ۷۸۵۲۵ رعمتاهال2۳60۵ 
لمعتاوتاهاه مطا بجمتامیی لهمزع1۲0۱۵ بط مط ص .07و۲۵ ,0.99 ,0.21 ,0.99 2208 0.99 ,0.16 ,0.99 0062168 ۷7۵۲۵ 
مطا 06۲ص افص 0صنام مه 16 ,افص .صمتاهلته۷۱ مه 4ج ۵۵00 ممتاهتطاتله مطا ج0 ومناله۷ ع1صما260800 ۷۵۲۵ 1106112 
6 0 1201 2 0 مب عصتصتامع0 مواق و1 ممتاهام‌وصهتام و۵72 امک .عمقهع0۲ع0 ۲۱۷۵۵۲۶ رعصقل ماقصتاه ۵۶ 1011۵86 
۰ ۱۷۱ ۲۵12010850 دمم عص مه ال مطلوهم10 ود ممتاهامممهاممه۱ه لفتاصمامم تاه رتم۳۷2 


٩۷۷۲۸۵۸۲ ۱۷۲۵061, ۷۷2۱۵ 6‏ ,1۳0۳/۵۹62 رومتتقوهعو ملقصتات رمعصقهه مفمصتلن 6۲۱۷۵0۸۵۵۵۰ 


